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Constructing regular self-similar solutions to the 3D Navier-Stokes equations 
originating at singular and arbitrary large initial data 

Zoran Grujic 1 



Abstract 

' Global-in-time smooth self-similar solutions to the 3D Navier-Stokes equa- 

tions are constructed emanating from homogeneous of degree —1 arbitrary large 
initial data belonging only to the closure of the test functions in Lf 
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1 Introduction 



A very simplistic description of the existing global-in-time results for the 3D Navier- 
Stokes equations is as follows: weak solutions are constructed from arbitrary large 
(mostly) L 2 -type initial data, while mild and strong solutions are generated by small 
initial data in a variety of functional spaces. It is not known if the weak solutions are 
regular and there are only partial regularity results - Caffarelli, Kohn and Nirenberg 
|CKNj proved that in a class of weak solutions satisfying localized energy inequal- 
ity, the so called suitable weak solutions, for every time T > the one- dimensional 
Hausdorff measure of the (possible) singular set in Q x (0, T) is (Q denotes the 
spatial domain). The standard constructions of weak solutions including the original 
Leray's construction [E] and the construction of the suitable weak solutions CKN 
start from the inital data belonging to L 2 on the whole spatial domain. Recently, 
Lemarie-Rieusset |L-Rlj (henceforth, the references will be to the book |L-R2j ) con- 
structed localized weak solutions started from initial data in the closure of the test 
functions in L 2 ocunif . 

Scaling-invariant solutions of any physical model are of special interest - heuristi- 
cally (and some times rigorously), they capture the qualitative behavior of the model. 
It turns out that there is a unique scaling invariance for the Navier-Stokes system: 
if u is a solution on M 3 x (0, oo) with initial data u$ then v x (x,t) = Xu(Xx, X 2 t) is 
a solution corresponding to the initial data Vq(x) = Xuq(Xx) for any A > 0. If the 
original initial data u are homogeneous of degree —1, i.e. if u (Xx) = X~ 1 u (x) for 
all x in M 3 and all A > 0, and the problem is considered in a uniqueness class for the 
Navier-Stokes equations, it follows readily that u(Xx, X 2 t) = X~ x u(x, t) for all x in M 3 , 
all t > and all A > 0, i.e. u is scaling-invariant or self-similar. It is worth noting 
that initial data generating (nontrivial) self-similar solutions are necessarily singular 
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- the homogeneity of degree —1 implies the existence of a singularity at the origin of 
the order of -rr. This prevents obtaining self-similar solutions in the 'usual' spaces, 
e.g., the Lebesgue spaces L p for p > 3 or the Sobolev spaces H s for s > 1/2. In fact, 
the j^j-type singularity is exactly a borderline singularity for the borderline spaces L 3 
and if 1 / 2 . The first rigorous construction of self-similar solutions is due to Giga and 
Miyakawa |GMj . They worked in the Morrey-type spaces of measures and obtained 
self-similar solutions for the vorticity equation starting from small initial data. More 
recently, Cannone [C] developed a systematic approach to constructing self-similar 
solutions in a variety of spaces (including the homogeneous Besov spaces) originating 
at small initial data. (See also [HKl IHMpI IHPl IFTl IF2j .) 

It is well-known there are no nontrivial self-similar solutions satisfying the Leray's 
energy inequality (see, e.g., |GMj ). Consequently, the standard constructions of weak 
solutions are not suitable for generating self-similar solutions. However, the weak 
solutions constructed in |L-R2j satisfy only the localized energy inequality and hence 
are not a priori incompatible with (nontrivial) self-similarity. 

In this paper we obtain existence of regular self-similar solutions emanating from 
arbitrary large initial data (homogeneous of degree —1) belonging only to Lf oc uni ^. 
To the best of the author's knowledge, this is the first time global-in-time existence 
of smooth solutions to the 3D Navier-Stokes equations is established without any 
smallness condition on the initial data. 

The proof starts off by showing that a modified local-in-time part of the con- 
struction of weak solutions in |L-R2j yields "partially self-similar" solutions on (0, T) . 
Partial self-similarity is then used to infer spatial regularity on (0, T) via Caffarelli- 
Kohn-Nirenberg estimate on the size of a singular set and local in space-time Serrin's 
L p x L q t -criterion. Spatial continuity in turn leads to full self-similarity on (0, T). This 
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interplay between self-similarity and regularity is exploited further until fully self- 
similar, C^C^-solutions on M 3 x (0, oo) are constructed. 

2 Preliminaries 

We consider the 3D Navier-Stokes equations (NSE) with the unit viscosity and the 
zero external force, 

u t - Au + (u ■ V)u + Vp = (1) 

V • u = 

where u denotes the velocity of the fluid and p the pressure. The spatial domain will 
be M 3 . 

In what follows, a projected form of the equation will be useful: 

u t - Au + PV- (u®u) =0 (2) 

V ■ u = 

(P denotes the Leray projector). Let us mention here that if a distributional solution 
u of (0) on M 3 x (0,T) is uniformly locally square-integrable, then there exists p in 
V'{R 3 x (0,T)) such that 

P V • (u <g> u) = V • (u ® u) + Vp 

in the sense of distributions, i.e. (u,p) is a distributional solution of (JTJ). 

Lemarie-Rieusset defined a localized version of weak solutions, the so called local 
Leray solutions. The definition is as follows |L-R2j . 
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Let Uq G Lj ocuni f with V • w = 0. A local Leray solution u corresponding to the 
initial data uq is a locally square-integrable distributional solution of (JIJ on M 3 x (0, T) 
with the following properties: 

i) uen t<T L°°((0,t),Ll niftloc ) 

ii) sup xoeR 3 n < s <t,\x-x \<i l v ® M ' 2 dxds <°° for ain < r 

iii) lim^o+ |w — Wo| 2 cfe = for any compact subset of M 3 

iv) m is a suitable weak solution in the sense of Caffarelli, Kohn and Nirenberg - 
in particular, u satisfies the localized energy inequality 

The following result can be found in |L-R2j - the idea of the proof is to use the 



localized energy inequality as the main tool in constructing the solutions originating 
at the L 2 OCjtmi j-initial data. Let E denote the closure of the test functions in the norm 
nf T 2 

Theorem 1. IL-RStf i) Let u G L 2 OCiUni f be divergence free. Then there exists a local 
Leray solution u onM. 3 x (0, T) for some T > 0. 



ii) If in addition uq G E, then u G H^tL 00 ((0, t), E) and linn,_+ + II tt — MoIIl 2 

loc,unif 



0. 



The construction is based on deriving some (uniform in e) local energy-type esti- 
mates for solutions of the family of mollified Navier- Stokes equations 

(u £ ) t - Au e + P V • {{u e * p e ) ® O = (3) 
V ■ u e = 

supplemented by u e (0, •) = u for all e > 0. A family of (standard) mollifiers {p e } e>0 
is defined as follows. Let p be a non-negative test function normalized such that 
J pdx = 1. Then p t is given by p e (x) = \p{-) for all x in R 3 . We will come back to 
this family of smooth approximating solutions in the following section. 



Remark 2. It was shown in |L-R2j that local-in-time local Leray solutions can be 
extended to global-in-time local Leray solutions; however, we will not need that result 
here. 

Next, we recall a definition of (space-time) singular /regular points of a weak solu- 
tion u in the sense of |CKNj . A point (x, t) in M 3 x (0, T) is singular if = oo 
for any space-time neighborhood D of (x,t). A point is regular if it is not singular, 
i.e. if there exists a space-time neighborhood D such that < oo. 

The main partial regularity result is given by the following theorem. 

Theorem 3. ICKM/ Let u be a suitable weak solution on IR 3 x (0, T). Then the 
one- dimensional Hausdorff measure of the possible singular set in M. 3 x (0, T) is 0. In 
particular, the singular set can not contain a smooth space-time curve. 

Remark 4. The original statement is somewhat stronger - it states that the one- 
dimensional parabolic Hausdorff measure of the singular set is 0. 

At first, it may be puzzling that local boundness is in this context identified as 
"regularity" . However, the following result due to Serrin [S] implies smoothness in x 
at any regular point (x,t). 

Theorem 5. JSJ/ Let u be a weak solution in some open space-time region D such that 
u e L p x L q t (D) for a pair of exponents (p, q) satisfying the Foias-Prodi-Serrin condition 
- + | < 1. Then u is of class C°° in x, and each derivative is bounded on compact 
subregions of D. 
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3 The construction 



Since our construction of global-in-time regular self-similar solutions will originate 
at the construction of local-in-time local Leray solutions, the first thing to check is 
whether the requirements on the initial data in Theorem Q are compatible with the 
desired self-similarity of solutions. As already pointed out in the introduction, any 
nontrivial self-similar solution must emanate from nontrivial initial data homogeneous 
of degree —1 and hence possessing a j^j-type singularity at the origin. This is not in 
conflict with being in Lf oc uni j, and decay at infinity, i.e. a possibilty of an approx- 
imation with the test functions, is also compatible with the homogeneity of degree 
— 1. A simple example of a family of homogeneous of degree —1 initial data satisfying 
the conditions of Theorem is the following: 



for a in E. Notice that this family generates arbitrary large initial data in Lf uni f. 

The second imminent problem in constructing any type of invariant weak solutions 
is that the weak solutions are not known to be unique, so although the inital data may 
be compatible with the invariance, it is a priori not clear if it is possible to construct 
any invariant solutions at all. One way out is to show that a particular construction of 
weak solutions will preserve the invariance - this was shown by Brandolese in the case 
of a special rotational invariance in the |CKJNj -construction of suitable weak solutions 
|Blj (as pointed out in [B2], the same argument works for any rotational invariance, 
and also for most of the standard constructions of weak solutions). The trick is to 
choose a test function in the definition of the smooth approximating solutions to be 
compatible with the invariance - in this case any radial test function will do. The 
consequence is the invariance of the approximating solutions which is then preserved 
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in the limit. The same trick does not seem to work for the self-similarity. It turns out 
that in order to construct self-similar approximating solutions, the test function would 
have to satisfy a certain degree of homogeneity which contradicts the very nature of a 
test function, i.e. having a compact support. What will happen here is that although 
the approximating solutions will not be self-similar, just enough self-similarity will 
be recovered in the limit to push the construction through. 

Theorem 6. Let uq G E be divergence free and homogeneous of degree —1. Then 
there exists a self-similar solution u of HJ) on M 3 x (0, oo) corresponding to the initial 
data uq with the following properties: 
i) limi^ + IK*) - u \\ L 2 =0 

iij u 6 C^C^° and hence a classical solution on M 3 x (0, oo) 
Hi) \\u(t)\\ L oo = 11^(1)11^ < oo for all t in (0, oo). 

Proof. Recall the construction of local-in-time local Leray solutions (jSJ). It is clear 
that will be scaling invariant if and only if the smoothed-out velocity 



exhibits the same scaling as the original velocity, i.e. if and only if w e (Xx, X 2 t) = 
X^w^x, t) for all A > 0. A simple computation yields w e (Xx,X 2 t) = A -1 f u t (x — 
z,t)X 3 p e (Xz) dz which is equal to X~ l w e {x, t) if and only if p e (Xz) = X~ 3 p e (z) for all 
A > 0, i.e. if and only if p e is homogeneous of degree —3. This is not possible since 
p € is just a rescaled test function and hence have a compact support. Consequently, 
an approximating solution u e will not be self-similar. 

However there is some symmetry. For 0<A<l,0<e<eo, define p^(x) = 
X 3 p e (Xx) and u^(x,t) = Xu e (Xx,X 2 t) for all x G M 3 and all t G (0,T). A straight- 
forward calculation reveals that u e solves (j2J) if and only if u^ solves the following 
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system: 



{u*) t -Au^ + PV- ((u x e * ® = (4) 



V • 11^ = 

supplemented by u x (0, •) = m (utilizing the homogeneity of u ). Since p^ ^ p t for 
any A ^ 1, the approximating solutions (for A ^ 1) and u e will differ as well. 
What saves the day is a simple observation that p^ = p± and hence = u±. Thus, 
{u t } is simply a subfamily of {u^}. However, we need to be somewhat careful here. 
Since a weak solution u in Theorem [T] is obtained through a sequence, say {u tn }, it 
is not possible to infer simultaneous convergence for continuum many lambda (and 
in particular, for A in (0, 1]). What we can do is to "preconditon" the sequence for a 
dense set. Let D 2 denote the set of all dyadics in (0, 1). Playing with indices in D 2 , 
we can construct (via a diagonalization procedure) a sequence {u En } with a property 
that {uzil} converges to a solution u for all A in D 2 U {1}. 

Recall that one of the convergences obtained in the proof of Theorem ^ is a 
strong convergence of ipu tn to <pu in LP ((0, T), L 2 (M 3 ) for any test function ip G 
T>(IR 3 x (0, T)), and any p < oo. Fix a ball B centered at the origin. Then, for any A 
in D 2 U {!}, a.e. t, 



u^ n (x,t)g(x) dx = / uzn_(x,t)g(x) dx — > / u(x,t)g(x)dx (5) 
Jb x Jb 

for all g in L 2 (B). On the other hand, it is easily seen that 

u^ n (x,t)g(x) dx = / Xu en (Xx, X 2 t)g{x) dx — > / Xu(Xx, X 2 t)g(x) dx (6) 
Jb Jb 

for all g in L 2 (B). Uniqueness of a limit coupled with (jSJ) and © yields that for any 
A in D 2 , a.e. t, 

[u(x, t) — Xu(Xx, X 2 t)~\ g(x) dx = 



B 



B 



for all g in L 2 (B). 

Now, fix A in D 2 and consider 

f(t)= / [u(x, t) — Xu(Xx, A 2 t)] g(x) dx 
Jb 

for t in (0, T). Since / is continuous (by the weak time continuity of u with values in 
L 2 (B)) and vanishes a.e., f is identically zero on (0, T). Next, for a fixed t in (0,T) 
consider 

g(X) = / £) — Aw(Ax, X 2 t)~\ g(x) dx 
Jb 

for A in (0, 1). This function is also continuous and vanishes on a dense set (D2) - 
hence vanishes for all A in (0, 1). Finally, since the ball B was arbitrary, 

u(x,t) = Xu(Xx, X 2 t) (7) 

for any A in (0, 1), any t in (0, T), and a.e. x in IR 3 . A priori, the a.e. x-set depends 
on both A and t. One can actually show more uniformity, but this will suffice for our 
construction. 

Notice that we can extend (J7|) to any A > provided that we restrict t to 
(0,min{T,£}). 

Next, we show that the singular set of u on R 3 x (0, T) is empty. In what follows, 
it will be convenient to define a family of parabolic space-time cylinders centered at 
a space-time point (x,t). For r > 0, define C x ^{r) by 

C Xjt {r) = {(y,r) 6l 3 x (0,oo)| \r - 1\ < K 2 , \(x,r) - (y,r)\ < r}. 

We argue by contradiction. Suppose there exists a singular point (x, t) of u in 
M 3 x (0, T) and consider a family of parabolic cylinders {C Xy t{r)} for < r < 
min{vr, \/T — t] = Rt,T- Since (x,t) is a singular point, (MIl^c^o-)) = 00 for 
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all r. Fix any A in ^0, t+ i R 2 j > an d then fix any r, < r < R t) r- For this choice 
of A, u(Xy, A 2 r) = A _1 w(2/, r) for any r in (t — |r 2 , t + |r 2 ), a.e. y in M 3 , and hence 
a.e (?/, r) in C^r). This implies that 



l M ll L °°(^, A 2 t (Ar)) =A 1 ||«|U»(C5 rit (r)) 



oo 



for any r, < r < -R^r- In other words, we constructed a family of parabolic 
neighborhoods of a point (Ax, A 2 t) on which u blows-up and so (Ax, A 2 t) is also a 
singular point. The argument can be repeated for any A in ^0, t+ ^ R 2 J an d hence 
we generated a smooth space-time curve of singular points passing through the point 
(x, t) contained in R 3 x (0, T). Since local Leray solutions are suitable, this contradicts 
Theorem El and thus the singular set must be empty. 

Now, it is standard to conclude infinite spatial regularity. Fix any space-time point 
(x, t). Since (x, t) is a regular point, there exists a (bounded) neighborhood D of (x, t) 
such that ||m||z,°°(d) < oo. This implies that u G L p x L q t (D) for any 1 < p, q, < oo, and 
in particular u satisfies the condition of Theorem It follows that u G C%?(D) and 
that each derivative is bounded on compact subregions of D. 

Let us note that the spatial continuity of u now implies that given A > 0, 

u(x,t) = Xu(Xx, X 2 t) (8) 

for all t in (0, min{T, J}) and all x G M 3 . 

Infinite time regularity in general does not follow from the emptiness of the sin- 
gular set. However, a restricted self-similarity will help. More precisely, we can write 

u(x, t) = —j=JJr ( —= J where Ut(- 

Since we know that u G C™, and it is clear that u G if and only if Ut G C°°, it 
follows that Ut is infinitely smooth. The infinite time regularity is now readily seen. 
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/I U {Jl '' l) f ° r aU ^ X,t) in R3 X ( °' T) - 



To summarize, u is in C^C^ (with the corresponding infinitely smooth pressure 
p) and hence a classical solutions on M 3 x (0,T) which is self-similar on (0,T), i.e. 
for a given (x, t) in M 3 x (0, T), 

u(Xx, X 2 t) = X~ x u{x,t) (9) 

for all A in ^0, yjjj (the corresponding pressure scales as p(Xx, X 2 t) = X~ 2 p(x,t)). 

Finally, we extend (u,p) to the whole space-time by self-similarity. It is clear 
that the extension will inherit infinite space and time regularity and will solve the 
equations (JJJ) in the classical sense (pointwise). 

One can now study various norms of u - we consider the L°°-norm here. Recall 
that u G H^t-L 00 ((0,t),E) which paired with the time continuity yields u(t) in E 
for any t in (0,T). In particular, u(j) is in E and hence it decays to as \x\ — > oo. 
Since we know u(f) is continuous, it follows < oo. Self-similarity implies 

IK*)IU- = \/| ||u (l) IIloo 7! 

for all t in (0, oo). This is informative at both ends. It gives a blow-up rate as t — » + 
compatible with the j^j-type singularity of the initial data, and also a generic decay 
rate as t — > oo. 

The convergence to the initial data in L 2 oc uni ^ follows from Theorem [TJ 

□ 

Remark 7. It is worth noting that most of the functional spaces utilized in the 
study of the 3D NSE on the whole space are embedded in L 2 oc uni j - u in L 2 oc gives the 
bilinear term V- (u®u) at least a distributional meaning, and uniformity corresponds 
to the (desirable) translational invariance of a norm. 
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Remark 8. It is clear that the proof of Theorem |U] implies infinite space and time 
smoothness of any self-similar suitable weak solution (regardless of its construction). 
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